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ABSTRACT 

Flows in which energy is transported predominantly as Poynting flux are 
thought to occur in pulsars, gamma-ray bursts and relativistic jets from compact 
objects. The fluctuating component of the magnetic field in such a flow can 
in principle be dissipated by magnetic reconnection, and used to accelerate the 
flow. We investigate how rapidly this transition can take place, by implementing 
into a global MHD model, that uses a thermodynamic description of the plasma, 
explicit, physically motivated prescriptions for the dissipation rate: a lower limit 
on this rate is given by limiting the maximum drift speed of the current carriers 
to that of light, an upper limit follows from demanding that the dissipation zone 
expand only subsonically in the comoving frame and a further prescription is 
obtained by assuming that the expansion speed is limited by the growth rate of 
the relativistic tearing mode. In each case, solutions are presented which give the 
Lorentz factor of a spherical wind containing a transverse, oscillating magnetic 
field component as a function of radius. In the case of the Crab pulsar, we find 
that the Poynting flux can be dissipated before the wind reaches the inner edge 
of the Nebula if the pulsar emits electron positron pairs at a rate N± > 10^° s~^, 
thus providing a possible solution to the "cr-problem" . 

Subject headings: pulsars: general - pulsars: Crab - MHD - plasmas 



1. Introduction 



Powerful outflows in which energy is transported predominantly as Poynting flux rather 
than kinetic energy arise in several astrophysical contexts. These include pulsars (Gunn 
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& Ostriker 1969; Rees & Gunn 1974), gamma-ray bursts (Usov 1992; Spruit 1999; Lyu- 
tikov & Blackman 2001), and the jets from galactic compact objects and active galactic 
nuclei (Lovelace et al 2002). An important generic property of such flows is that they are 
highly relativistic. Force-free solutions, valid in the inner regions, indicate that radial flows 
accelerate quickly (Buckley 1977; Contopoulos & Kazanas 2002) until they cross the fast 
magnetosonic point, at which their bulk Lorentz factor 7 satisfies 7 = \/a. The parameter 
(T, which is the ratio of Poynting flux to particle-born energy flux is by deflnition large under 
these conditions. Eventually, such flows release their energy into particles and/or radiation. 
In principle, this can happen if the flow diverges sufficiently rapidly and is thus forced to 
accelerate. However, this occurs only very slowly for ideal relativistic MHD flows which are 
initially radial (Beskin ct al 1998; Chiueh, Li & Bcgclman 1998; Bogovalov & Tsinganos 1999; 
Bogovalov 2001; Lyubarsky & Eichler 2001; Lyubarsky 2002a). In the case of non-steady 
winds, such as those driven by an obliquely rotating magnetic dipole, the current required 
to support the oscillating component of the magnetic field decreases more slowly with radius 
than the density of available charge carriers (Usov 1975; Melatos & Melrose 1996). Conse- 
quently, the ideal MHD approximation must break down at some point in the outflow where 
the implied minimum drift speed reaches that of light. However, internal dissipative effects 
may intervene before this stage is reached, reducing the wave amplitude and decreasing the 
required current. Shock fronts, which provide efficient dissipation for low a flows, are rather 
slow to dissipate the energy of a Poynting-flux dominated outflow, although they may be 
effective if the particle density is high (Lyubarsky 2002b). An alternative is dissipation 
by annihilation or reconnection of the oscillating part of the magnetic field (Coroniti 1990; 
Michel 1994), and this is the process we investigate in detail in this paper. The effectiveness 
of such a dissipation mechanism is an important ingredient in models of pulsed high energy 
radiation from pulsar winds (Kirk et al 2002; Skjaeraasen et al. 2002; Skjaeraasen & Kirk 
2001) and the after-glows of gamma-ray bursts (Drenkhahn & Spruit 2002). 

Although now invoked in models of AGN jets and gamma-ray bursts (Thompson 1994; 
Schopper et al 1998; Heinz & Begelman 2000; Spruit et al 2001), the original motivation for 
considering reconnection in Poynting-flux dominated outflows was the so-called "cr-problem" . 
This term is generally used to refer to a puzzle posed by the Crab Nebula and the pulsar 
wind which powers it: according to estimates of the number of pairs emitted by the pulsar 
(Arons 1983; Hibschman & Arons 2001b; Gallant et al 2002), the energy flux carried by the 
wind whilst close to the star is dominated by Poynting flux (Michel 1982). However, at the 
inner edge of the Crab Nebula, where the wind passes through a termination shock, it must 
have a very low magnetization, a ~ 10^'^, if it is to match the observed expansion speed of 
the nebula (Rees & Gunn 1974; Kennel & Coroniti 1984; Emmering & Chevaher 1987) and 
the axial ratio of its optical isophotes (Begelman & Li 1992). Thus, unless the structure of 
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the termination shock is very different to that of a hydrodynamic shock, dissipation of the 
magnetic field in the fiow must occur before the fiow reaches this radius. 

The wind of the Crab pulsar is driven by a rotating neutron star in which a non- 
axisymmetric magnetic field is embedded — often assumed to be an oblique magnetic dipole. 
Close to the equatorial plane, the magnetic field in the wind forms toroidal "stripes" of 
opposite polarity, separated by a current sheet (Coroniti 1990; Bogovalov 1999), as illustrated 
in Fig. 1. In principle, a dissipative process operating on the scale of the wavelength of the 
striped pattern is capable of annihilating the oscillating component, which, in the equatorial 
plane, corresponds to the entire magnetic field. But the rate at which dissipation proceeds is 
uncertain. Using a particular prescription in which the speed of the charge carriers attains 
the maximum value of c, Coroniti (1990) and Michel (1994) found for the case of the Crab 
pulsar that dissipation should be completed well inside the termination shock, which occurs 
at a radius of roughly lO^rL (where tl = c/Q, is the radius of the hght-cylinder and Q, the 
angular velocity of the star). However, an important effect of dissipation is that it causes the 
fiow to accelerate. As magnetic field energy is converted into heat, the plasma performs work 
on the expanding wind and converts the heat into kinetic energy. Because of the relativistic 
dilation of time, this effect, which was not included in the treatments of Coroniti (1990) and 
Michel (1994), leads to an apparent decrease in the rate of dissipation in the rest frame of 
the star [Lyubarsky & Kirk (2001) henceforth LK]. As a result, this rate of dissipation does 
not solve the cr-problem. 

Physically, this approach to dissipation is based on a one-dimensional picture in that 
it corresponds to magnetic field annihilation in a current sheet that is uniform in the plane 
perpendicular to the radius vector. In reality, we might expect the sheet to break up into a 
large number of regions where dissipation proceeds via magnetic reconnection in fiow patterns 
which are two or three dimensional in nature (see, for example Priest and Forbes 2000). In 
this case, the global effect on the fiow can be found by (i) maintaining the assumption that 
dissipation is confined to a layer around the field reversal, but allowing the thickness of this 
layer to be determined by an estimate of the average dissipation rate based on the physics of 
reconnection and (ii) assigning to the layer an average energy density and pressure, related 
to each other by an effective equation of state. The prescription used by Coroniti (1990), 
Michel (1994) and Lyubarsky & Kirk (2001) is equivalent to setting the thickness of the 
dissipation layer equal to the gyro radius of the current carriers and using the equation of 
state of an ideal relativistic gas. This gives a reasonable lower limit, but other prescriptions 
have been proposed (Lyubarsky 1996; Drenkhahn 2002; Lyutikov & Uzdensky 2002). In this 
paper we argue that an upper limit to the reconnection rate can be implemented by limiting 
the expansion speed of the layer to that of sound in an unmagnetized gas. We present the 
global properties of reconnecting, spherically symmetric, relativistic winds, including the 
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Fig. 1. — Top: The Parker spiral pattern in the equatorial plane. The arrows indicate the 
magnetic field direction of an obliquely rotating magnetic dipole, and the solid line gives the 
position of the current sheet. Bottom: the poloidal structure of the current sheet (thick, 
solid line). The arrows denote the spin axis (vertical) and magnetic dipole axis (inclined), 
and the thin line bounds the equatorial zone of the wind, in which magnetic reconnection is 
possible. 
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dependence of the Lorentz factor on radius, for each prescription. Our findings apply to any 
spherical relativistic flow that contains a toroidal field with radial fluctuations, but we phrase 
the discussion in terms of a pulsar wind, since this fixes the wavelength of the fluctuation 
as well as its initial amplitude. In the particular case of the Crab pulsar, we show that the 
a-problcm can be resolved if the initial value (i.e, the value at the fast magnetosonic point) 
of this parameter is at most 10^, or, equivalently, if the /^-parameter [see Eq. (23)] introduced 
by Michel (1969) is less than a few times 10^. 

The paper is organized as follows: in Sect. 2 we review the equations derived by 
Lyubarsky & Kirk (2001) that describe, in the short wavelength limit, spherically symmetric 
MHD winds containing an oscillating component of the magnetic field which can undergo 
damping by dissipative effects. In Sect. 3 we discuss the modelling of the dissipative layer 
and describe how it can be incorporated into the MHD description. Section 4 contains a 
discussion of the various prescriptions for dissipation. In section 5, asymptotic solutions, 
valid for relativistic plasma temperatures, are found that give good approximations to the 
global properties of the flows, such as the radius at which dissipation is complete. These are 
compared to numerical solutions, which indicate the sensitivity to variations in the initial 
conditions, delineate regions in which the dissipation prescriptions become unphysical and 
quantify departures from the asymptotic solutions when the plasma temperature is nonrel- 
ativistic. We compute the effective diffusivity implied by our results in Sect. 6 and compare 
it to the Bohm diffusivity. The implications of our results, in particular for the Crab pulsar 
and its nebula, are discussed in Sect. 7 and a summary of our conclusions is presented in 
Sect. 8. 



2. MHD equations in the short wavelength approximation 

In this paper we investigate the evolution of a radial MHD wind containing a purely 
toroidal magnetic field. This idealisation is thought to be a reasonably good approximation 
in the case of a pulsar for radii far outside the light cylinder, r ^ tl [see, for example Coroniti 
(1990)]. However, alternative pictures in which the radial component of the magnetic field 
plays a role are possible (Kuijpers 2001). Following Lyubarsky & Kirk (2001) we consider the 
slow evolution of a periodic entropy wave, containing two reversals of the magnetic direction 
per oscillation, as it propagates outwards in the spherically symmetric wind. The flow pattern 
is assumed to be steady in the sense that at fixed radius the amplitude of the wave is constant 
in time, so that the flow variables are exactly periodic. The radial evolution is governed 
by two factors: the gradual annihilation of the oscillating component of the magnetic field 
embedded in the wave according to one of the microphysical prescriptions discussed in Sect. 4 
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and the expansion and acceleration of the plasma as a result of the spherical geometry and 
the outwardly directed force arising from the pressure gradient. The evolution of the wave 
is assumed to be slow — significant changes are permitted only over radial scales large 
compared to the wavelength (= 27rrLf /c, where v is the bulk radial velocity). The equations 
describing the evolution are obtained from a perturbation expansion in the small parameter 
r-i^/r and involves averaging over a wave period the equations of continuity, energy and 
entropy. The procedure is described in detail in LK. The zeroth order equations give the 
familiar conditions for the entropy wave mode: pressure equilibrium between plasma and 
magnetic field and zero bulk speed in the wave frame. Imposing non-secularity conditions 
on the first order equations then gives expressions governing the slow evolution of the zeroth 
order quantities. Prom the continuity equation one obtains: 



|-(rV(n')) = 0, 



and from the energy equation: 



471-^2 



= 0. (2) 



The entropy equation includes ohmic dissipation to first order in ri^/r, and can be written 
in terms of the magnetic field, using Maxwell's equations: 

{e' + p') d , 2 N 

In these expressions, gaussian units have been employed, the magnetic field and fiuid velocity 

— * 

— * — * 

in the laboratory frame (in spherical polar coordinates) are B — B{r,t)(l) and v — v{r,t)f 
and the notation (•) denotes an average over a wave period. The proper particle density is 
denoted by n', the proper energy density e' is the component T^'" of the stress-energy tensor 
measured in the rest frame and the pressure p' refers to the radial component T^'^ in the 
same frame. 



3. The current sheet 

In order to perform the averaging in Eqs. (1), (2) and (3) an explicit model of the 
entropy wave and the associated current is needed. Lyubarsky & Kirk (2001) used a fluid 
model in which a hot, unmagnetized plasma separates regions of opposing magnetic field 
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containing cold plasma. The transition at the interface of the hot and cold regions was 
assumed to be sharp. In this paper, we use instead an exact stationary solution of the full 
Vlasov/Maxwcll kinetic equations that contains a magnetic field reversal. This solution also 
forms the basis of the discussion of the microphysics presented in Sect. 4. When an average 
over one wavelength is performed, the equations obtained are essentially the same as those 
used by Lyubarsky & Kirk (2001), with two minor differences. The first is connected with 
the presence of a cold plasma component within the hot parts of the sheet and is noted below 
where it appears. The second arises when the drift speed of the charge carriers is relativistic 
and it becomes necessary to compute the effective ratio of specific heats of the plasma as a 
function of position according to Eq. (A23). We expect this effect to be unimportant and do 
not investigate it in detail. The use of an equilibrium distribution presupposes an effective 
collision time tcoU that is short compared to the dynamical time tdyn- On the other hand, 
each dissipation prescription has an associated value of tcou- In Sect. 6 we show that indeed 
^cou < ^dyn in each case. 

The exact stationary solution we use was originally found by Harris (1962) for the 
nonrelativistic case and subsequently generalized to apply also to the case of relativistic 
particle temperatures and drift speeds by Hoh (1968). It contains three parameters which 
can be chosen to be the magnetic field B'q far from the current sheet, the number density of 
particles of each charge at the center of the sheet and the sheet width a, each measured 
in a frame of reference in which the sheet is stationary and the electric field vanishes. One 
can then define a temperature T (in energy units) which is independent of position: 



The Debye length at the center of the sheet 



is, to within a factor of the order of unity, equal to the gyro radius of a particle of energy 
T in the magnetic field B^, provided T ^ mc^. The sheet is neutral and the electric 
field vanishes everywhere. The current is carried by oppositely directed drifts of oppositely 
charged particles. Specialising to the case in which only positrons and electrons are present, 
the speed ±c/3 with which each component drifts is 



and is independent of position in the sheet. The current sheet in this solution is planar, but 
can be used to describe each of the sequence of nested quasi-spherical sheets contained in the 




(5) 




(6) 
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radial wind, provided the sheet thickness remains small compared to the wavelength. The 

frame of reference used to portray the sheet solution is one which moves radially outwards 
with the bulk velocity of the wind. A derivation of the internal structure is presented in 
Appendix A. For a single sheet located at radius r at t = 0, the magnetic field B'[r) and 
the density of each species of charge carrier are 



B'{r,t) = B'o{r) tanh 
n'^{r,t) = A^^(r)sech^ 



7ft 
a(r) 

•yvt 



(7) 
(8) 



a(r) 

An additional spatially uniform component of cold, neutral, current-free plasma of arbitrary 
density can be added to this solution without changing the above results. The solution is 
also unaffected by the addition of a constant magnetic field of arbitrary strength in the 9 
direction. 

One wavelength of the entropy wave contains two current sheets of opposite polarity. 
Defining the fraction A of the wave occupied by hot plasma as 

A = (9) 

and averaging over one wave period, one finds for the particle number density n' {—2n'_^): 

{n') = 2iV4A + <, (10) 

where n'^ is the density of the uniform cold background plasma. The average of the pressure 
p' follows from the stress-energy tensor of the sheet particles (Eq. A21) 

{p') = 2N'^AT = I^A, (11) 

and this can be related to the average energy density by defining the effective ratio of specific 
heats 7 according to 

p' = (7 - 1) (e' - n'mc^) , (12) 

[see Eq. (A23)]. In principle, 7 is capable of accounting not only for potential anisotropics 
in the gas {p' refers to the r-r component of the stress-energy tensor — see Appendix A) 
but also for other non-equilibrium particle distributions, e.g., those of power-law type. 

The magnetic field as given by Eq. (7) leads to 

B^\ _ /B^' 
^ 4717^ / \ 47r 
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d 

B— irvB) 
or 



/2 



An 

..2 D/2 



A), 



d 



r^o (1 - A) ^ (rv) + 



rv d ' 



I'B': (1 - A) 



(13) 



(14) 



Using these expressions, the (integrated) continuity and energy equations and the en- 
tropy equation (1, 2 and 3) become: 



n/2 

2^ I 2 2 -'-'O 

7mc $ + r 7 v 

Stt 



2 + 



7 



7 



7 



2iV± 
L 

7 - 1 y J ' 



(15) 
(16) 



9 



7U— 



A 5'^ 



9r \7 — 1 Stt 



V d 
7 (9r 



7^(1 



A) 



Stt 



7^^,^ (1 - A) a , , 
+ ^ : -— {rv) 



Airr 



dr 



0, 



(17) 



where N± is the rate at which pairs leave the star (corresponding to 2N± particles per 
second) and L is the luminosity in the wind, which is assumed to occupy a solid angle Q^. 
The constants $ and ^ represent the particle flux and the total energy flux. 

These equations are complemented by the condition that far from the center of the 
current sheets the magnetic flux remain frozen into the plasma: 



rn' 



constant. 



(18) 



We now define the fraction ^ of particles carried in the hot plasma: 

27V4A7wr^ 



(19) 



and a dimensionless quantity p that is related to the magnetic pressure in the cold part of 
the wind: 



P = 



r^-i^vB'^^^ 
47r^ 



(20) 
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The continuity equation then simply relates the density of cold plasma n'^to 



and the energy equation becomes 



7 



+ 



1 + 



7 



P 



// L 2(7-1) 

where the constant /i is the total energy carried by the wind per unit rest mass (xc^): 

* L 



(21) 



(22) 



(23) 



mc2$ 27V±mc2 ' 

as introduced by Michel (1969). Physically, (1 — A)p is the fraction of the wind luminosity 
carried by magnetic field, 7//^ is the fraction carried by the particle rest mass and the 
remainder, p7A/[2(7 — 1)], can be attributed to the enthalpy flux. The constant ji is related 
to the more commonly used a (Kennel & Coroniti 1984; Kirk & Duffy 1999), defined as the 
ratio of Poynting flux to particle-born energy flux by 

-1 



a = (l-A)p 



7 



7-1 



(24) 



where a is, in general a function of r, but is constant in cold (A = 0), spherically symmetric 
flows of constant speed, in which case a — ^ip/^ ~ /x/7, for ^ 1. 



The entropy equation is 

'2-7 



7 



dA 



+ 



1 + 



+ 



1 / d In r 

c2 /2-7 



7 



7-1 



7-1 
d ln7 
dlnr 



dlnp 
dlnr 

+ 2A = 0, 



(25) 



in which, for simplicity, we have assumed 7 is constant, and the condition of flux freezing in 
the cold part of the flow is 



pv 



constant. 



(26) 



Thus the conservation equations and the entropy equation yield three equations (22), 
(25) and (26) in the four unknowns ^, 7, p and A. ^ 



-'^The current sheet model adopted by Lyubarsky & Kirk (2001) does not allow cold and hot components 
of the particle distribution to coexist at the same position. This leads to a slightly different expression for 
the density of cold particles: ^ + ^vr'^n'^{l — A)/$ = 1, instead of Eq. (21), and, consequently for the flux 
freezing condition: p (1 — A)^ v/{l — ^)^ = constant, instead of Eq. (26). 
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To close this set, it is necessary to specify the rate at which dissipation occurs using one 
of the prescriptions discussed below. An additional physical constraint on these can be found 
by requiring that the ohmic dissipation term in the entropy equation be positive. Performing 
the perturbation expansion as indicated above and described in Appendix B, this gives the 
condition 

> 0. (27) 



dlnr 



P 



1/7 



2\(7-l)/7 



4. Dissipation prescriptions 

In the following, three different prescriptions for the effective, globally averaged rate of 
dissipation are introduced. Each prescription has a specific, physical motivation and corre- 
sponds in a resistive MHD analogy to a certain value for the effective magnetic diffusivity, 
as discussed in Sect. 6. 



4.1. Slow dissipation 

The current sheet described above can be realized only if the drift speed required by 
Eq. (6) is subluminal. This gives an absolute lower hmit on the thickness of the sheet, or 
alternatively, on the fraction of particles it contains. The Debye length can be written as: 

Ad /.r fpvy'^^ ^28) 



where we have defined the dimensionless fiow luminosity 



L = I'S 1 . (29) 



Then the condition for a subluminal drift speed, (3 < 1, becomes 

e>e.,. = ^ff)"". (30) 

rLl \LcJ 

The treatments of both Coroniti (1990) and Michel (1994) amount to setting C, = C,mhi- 
Lyubarsky & Kirk (2001) also used this condition. A current sheet with particle distributions 
that are in equilibrium does not, of course, dissipate magnetic energy. However, as the 
sheet moves outwards with the fiow, the parameters describing it change and it is this slow 
evolution through a sequence of quasi-equilibria which implies dissipation. Eq. (30) leads to 
a lower limit on the dissipation rate, since a larger value of ^ implies a thicker sheet (larger 
A) and smaller fraction of the total energy fiux carried by the fields. 
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4.2. Fast dissipation 

The maximum dissipation rate is more difficult to estimate. It follows from the max- 
imum speed at which the current sheet can expand into the surrounding plasma and still 
remain in the neighbourhood of the Harris equilibrium. Denoting this speed by cj3c (> 0), 
the associated rate of evolution of the sheet thickness with radius, can be written: 

dA 2l3^rc^ 



dlnr 7rrL7^ {v"^ — c^Pc) ' 

(see Appendix C). 



(31) 



In the outer parts of the sheet, the balance between magnetic and particle pressure can 
no longer be maintained if the bulk velocity exceeds the local speed of the fast magnetosonic 
mode c/5f. This suggests the condition 



/3c < A = 



5'V47r 



1 + ^(1-0 
HP 



1-1/2 

(32) 



Drenkhahn (2002) has generalized this prescription, introducing a parameter e = /3c/A ^-iid 
arguing that it is typically of the order of unity. 

— * 

However, in the absence of a component of the magnetic field in the 6 direction, the 
Harris equilibrium has at its center a neutral layer. This can only be maintained close to 
equilibrium if the expansion speed is small compared to the sound speed in the unmagnetized 
gas. Were this not the case, the layer would lose its internal pressure support and be crushed 
by the external magnetic pressure. In the case of a relativistic temperature (T ^ mc^) this 
leads to the condition 

A < (33) 

In the outer parts of the flow, where the plasma temperature decreases to below the electron 
rest mass, condition (33) is inapplicable and for T <C mc^, the sound speed is \J^T j {m(?). 
A full treatment of the dynamics in this regime requires a local evaluation of the equation of 
state, i.e., of the effective ratio of specific heats 7. In this paper, we avoid such complications 
by taking 7 = 4/3 everywhere in the current sheet, which is appropriate for any isotropic 
gas in which the particles are relativistic [E ~ cp). Nevertheless, we can take approximate 
account of a more realistic sound speed by setting 
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In flows in which the Poynting flux dominates, this restriction leads to a much more stringent 
upper limit on the dissipation rate than does Eq. (32), since /3f ~ 1, for a ^ 1, and the 
denominator on the right-hand side of Eq. (31) becomes small. The limit can be relaxed 



the sheet. In this case the local speed of the fast mode would be more appropriate than the 
speed of sound in an unmagnetized plasma. However, this makes only a minor difference, 
provided the plasma at the center of the sheet is not too far from equipartition. 



The Harris equilibrium is well-known to be unstable to the excitation of a variety of 
waves in the non-relativistic case — see, for example, Daughton (1999); Biskamp (2000). 
It is thought that the current sheet first filaments because of the growth of the tearing 
mode, producing magnetic islands separated by neutral lines. These then lead to magnetic 
reconnection and its associated dissipation (Melrose 1986). This scenario has largely been 
developed in the context of the geomagnetic tail of the Earth. Particle-in-cell simulations 
appear to confirm its vahdity (Zhu & Winglee 1996; Pritchett et al. 1996), although the 
3-dimensional picture is more comphcated (Biichner & Kuska 1999; Hesse et al. 2001) 

In the relativistic case the linear growth rate of the tearing instabihty has been calculated 
by Zelenyi & Krasnosel'skikh (1979). For T » mc^ they find for the most unstable mode, 
which has a wavelength approximately equal to the sheet thickness. 



Thus, for a given number of particles in the sheet, and fixed external pressure (implying ^, 
p and /3 constant), the growth rate is faster for thinner sheets. Lyubarsky (1996), in his 
model of the high energy emission from the Crab pulsar, suggested that the speed /3c of 
expansion of the sheet could be self-regulating being governed by the fastest growing mode 
of the relativistic coUisionless tearing mode instability: 



slightly if an additional component of the magnetic field (in the d direction) is present in 



4.3. Dissipation limited by the teeiring mode 



7t - P'Na, 



(35) 



/3c = a7t/c. 



(36) 



or, in terms of the fiow variables. 




dlnr 




(37) 
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5. Solutions of the flow equations 

In this section we estimate both the length scale on which dissipation converts an initially 
Poynting-dominated flow into one in rough equipartition, as well as the rate at which the 
flow accelerates, i.e., 7(r) during this process. To simplify the results, we specialize to the 
case of a relativistic gas: 7 = 4/3. 

At an initial radius, assume the flow contains a thin sheet A ^ 1 containing relatively 
few particles ^ <^ 1, is Poynting flux dominated with 7//^ -C 1 and already super magne- 
tosonic 7 > /x^/^. As the radius increases, the flow accelerates, but, as long as the above 
ordering is maintained, we may set p ~ 1 and v c. In this range, the flux freezing equation 
(26) reduces to the approximate relation 

- - 2^ + A = constant, (38) 



and the entropy equation is 



- , 2A + ^ ) -— In7 + 2A = 0. (39) 
dmr \ /x/dmr 

When these two equations are complemented by one of the dissipation prescriptions 
[using, for fast dissipation, Eq. (33)], it is straightforward to obtain asymptotic solutions for 
the flow for large radius. These take the form 

A oc r'', (40) 

e ^ ^. («) 

1 = (42) 

Id q 

[or, in the case of the LK sheet, ^ = A(l+g)/g instead of Eq. (41)]. Shortly before dissipation 
is complete, at the point where p decreases, A ^ ,^ ^ 1 and ^ fi, the solutions lose their 
validity. For each dissipation prescription, the asymptotic solutions yield reasonably accurate 
estimates of the radius at which this occurs, as well as the dependence of the flow variables on 
radius up to this point. An estimate of the smallest radius at which the asymptotic solutions 
are valid is obtained by locating the magnetosonic point. These results are summarized in 
Table 1. The radius at which dissipation is complete can be estimated by inserting into the 
asymptotic solution either A = l,^ = l, or7 = yU and the numerical integrations indicate 
that the latter is somewhat more accurate. Table 1 gives this radius with a numerically 
estimated coefficient to order of magnitude accuracy. For slow dissipation this asymptotic 



- 15 - 




. I . I . I . I . I . I . I 

2 4 6 8 10 12 14 

log(r/rJ 




log (r/rj 




I I I I I I . I . I I I I 

2 4 6 8 10 12 14 

log (r/rj 




log (r/rJ 




I I I 



2 4 6 8 10 12 14 

log(r/rJ 



Fig. 2. — The flow variables as functions of radius for slow dissipation and parameters 
L = 1.5 X 10^^, f2w = 47r, yU = 2 X 10^, thought appropriate for the Crab pulsar. Displayed 
are the average entropy [see Eq. (27)], the bulk Lorentz factor, the fraction of a wavelength 
occupied by the hot sheet A, the temperature and the Mach number of the expansion speed 
of the sheet. The top right panel shows the ratio ^/A (which gives the ratio between the 
density of hot particles in the sheet, and the wavelength-averaged density of both hot and cold 
particles). Three numerical solutions are shown as solid lines for different initial conditions 
on A. The asymptotic solutions for 7 and A are shown as dashed lines. 
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Fig. 3. — The flow variables as functions of radius for tearing mode limited dissipation 
(parameters as in Fig. 2). Three solutions are shown as solid lines for initial conditions 
A = 10~^ and 10~^ starting at r = lOrL, and for A = 10~^ starting at r = lOOOrL. The 
asymptotic solution is shown dashed line. 
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Fig. 4. — The flow variables as functions of radius for fast dissipation according to Eqs. (31) 
and (33), for the same parameters and initial conditions as in Fig. 3. The solutions with 
initial values of A = 10~^ and A = 10"'' coincide for almost all radii. 
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solution was originally found by Lyubarsky & Kirk (2001); for fast dissipation, the scaling 
with radius agrees with that given by Drenkhahn (2002). 

Figure 2 compares the results of a numerical integration of the full equation set, Eqs. (22), 
(25) and (26), together with the prescription for slow dissipation in which ^ is required to 
equal ^min, defined in Eq. (30). The parameters chosen in this and the following three fig- 
ures arc those implied by the discussion of the Crab pulsar in Sect. 7: L — 1.5 x 10^^ and 
// = 2 X 10^. Six quantities are plotted against radius: 

1. the quantity defined in Eq. (27) which is proportional to the entropy averaged over a 
wavelength 

2. the ratio ^/ A, that is the ratio of the density of hot particles to the wavelength-averaged 
particle density 

3. the Lorentz factor 7 (in units of its maximum value /i) 

4. the fraction A of a wavelength occupied by hot plasma 

5. the temperature 

6. the absolute value of the Mach number (with respect to the fast magnetosonic speed 
in the cold part of the wind) with which the sheet expands. 

Solutions are shown for three different sets of initial conditions with A = 10~^, 10~^ and 
10~^, all starting at the point r = lOrL with a slightly supermagnetosonic Lorentz factor 7 = 
1.1/i^/^ 33. All converge rapidly to the asymptotic solution given in Table 1 and Eqs. (40- 
42), which, in the case of 7 and A, is plotted as a dashed line. Deviations from this solution 
are small and appear only close to the maximum radius, where dissipation is complete. Using 
this prescription for dissipation, it is possible to choose an arbitrarily small initial value of 
A — the results do not differ appreciably from the case A = 10^^. However, these solutions 
adjust rapidly in the initial phase and cannot strictly be described by the perturbation 
approach underlying our method, which assumes all variations to be on a length scale equal 
to or larger than the local radius. For the larger initial values of A, the region occupied by 
hot plasma shrinks as the flow moves outwards. This is accompanied by a decrease in the 
averaged entropy, indicating that the prescription is unphysical. Thus, we cannot be sure 
that initial states with such a large value of A will approach the asymptotic solution, since 
a different dissipation prescription with increased entropy generation is required at least 
initially. Different values of the starting point in radius yield qualitatively similar solutions. 
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Fig. 5. — The flow variables as functions of radius for fast dissipation using a realistic 
estimate of the sound speed according to Eq. (34), and for the same parameters and initial 
conditions as in Fig. 4. The solutions with initial conditions A = 10~^ and A = 10~^ coincide 
for almost all radii. 
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Figure 3 shows results for the case of tearing mode hmited dissipation, Eq. (37). Three 
different sets of initial conditions were chosen: A = 10~^ and 10~^ starting at r = lOrL and 
A = 10~^ starting at r = lOOOrL- Because this prescription takes the form of a first order 
differential equation, an additional initial condition on C, is needed — we choose ^ = A. Once 
again, the asymptotic solution provides a reasonable approximation at large radius. In this 
case, however, the entropy increases in all solutions, despite the fact that the region occupied 
by hot plasma shrinks slightly. Thus, this prescription provides a physically consistent 
description of the approach of the flow to the asymptotic solution even for quite large initial 
values of A. However, the initial density contrast ^/A between the hot and cold parts of the 
flow has to be chosen with with some care in order to avoid rapidly varying transients. At 
r > lO^rL, the plasma temperature is nonrelativistic, T < rru? and the expression used for 
the growth rate of the tearing mode instability loses validity. 

The case of fast dissipation is shown in Fig. 4, where we have used the relativistic sound 
speed (33) as an upper limit on the expansion rate of the sheet. Here again, the dissipation 
prescription takes the form of an ordinary differential equation, so that initial conditions 
on both A and ^ are required. These were chosen to be the same as those of Fig. 3. The 
solutions again converge to the asymptotic expressions, and the entropy rises throughout. 
However, the initial transients fluctuate more strongly here than in the case of the tearing 
mode limited prescription. The result of the more realistic estimate of the sound speed given 
by Eq. (34) is shown in Fig. 5, for the same parameters and initial conditions as used in 
Fig. 4. Differences between these two sets of solutions appear only in the outer parts of 
the flow, where the temperature of the plasma becomes nonrelativistic. The Mach number 
drops sharply at this point, and the dissipation rate is reduced. Independent of the initial 
conditions, the solutions terminate at r 2 x lO^rL. No asymptotic solutions are available 
to describe this outer phase of the flow. Nevertheless, the estimates of the radius at which 
dissipation is complete obtained using the asymptotic solutions for fast reconnection given 
in Table 1 remain accurate to within roughly an order of magnitude. 



6. Compcirison with dissipation by Bohm diffusion 

For each dissipation prescription, the effective relaxation time of the plasma can be 
estimated by appealing to an MHD picture: the magnetic diffusivity r] is related to the sheet 
thickness and the dynamical evolution time in the comoving frame tdyn by 

ri ^ f- (43) 
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and is given in terms of the effective collision time icou by 

ri - ^, (44) 

where we have used Ad instead of the Larmor radius, as appropriate in a relativistic plasma 
(in the current model the two differ by a constant factor of the order of unity). Combining 
these two equations gives 

^coU ~ ^dyn (45) 

and we see that the collision time is shorter than the dynamical time provided a > Ad, or 
equivalently, ^ > ^min- 

As can be seen from Eq. (40), the asymptotic solutions presented in Sect. 5 correspond 
to a dynamical timescale for the current sheet in the comoving frame given, to within a 
factor of the order of unity, by 

idyn = — ■ (46) 

Thus, it is straightforward to compare the effective collision time implied by the dissipation 
prescription ^cou with the relaxation time ~ Ad/c corresponding to "Bohm" diffusivity in 
a relativistic plasma. Using Eqs. (45) and (28), we have 

^cou ^ /Ad 

V « y V^lA72 

C^^rlA [lcJ ■ ^ ^ 

At the outer radius, this ratio takes on its smallest value for slow dissipation and for dissi- 
pation limited by the tearing mode: 



^ [ L^' ^/ /J, j (tearing mode) 



and its largest value for fast dissipation: 



-2 

(fast). (49) 
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7. Discussion 

For the acceleration of the flow to proceed as we have described, the actual dissipation 
rate must lie above that of 'slow' dissipation and below that of 'fast' dissipation. To a reason- 
able approximation, these conditions can be imposed simply by inspecting the appropriate 
values of rmax and rmin, as given in Table 1. It is immediately clear that no prescription can 
fulfil these requirements unless // is less than a certain maximum value. Adopting /3c = 1/ -\/3, 
we find 

H < SL^/^ . (50) 

Physically, restriction (50) arises because at higher values of /j, the sheet would have to expand 
supersonically in order to accumulate enough particles to maintain the required current. 
This restriction was pointed out by LK, who speculated that, when it is violated, dissipation 
ceases before the available magnetic energy is exhausted. Whether this is the ultimate fate 
of such a wind is not completely clear. It is certainly true that our treatment using the 
Harris sheet or a model with hot and cold fluids fails in this regime, and it seems probable 
that no equilibrium can then exist in which the electric field vanishes in the comoving frame. 
Perhaps at this point it is essential to allow the conversion of the entropy wave into an 
electromagnetic mode, along the lines suggested by Melatos (2002). 

If Eq. (50) is fulfilled, Eqs. (48) and (49) show that our solutions imply slower dissipation 
than the Bohm value at all radii for the 'slow' and 'tearing-mode' cases, but faster-than- 
Bohm dissipation in the 'fast' case. In standard treatments of reconnection, this value is 
often taken as an upper limit for the diffusivity (Somov & Titov 1985; Priest & Forbes 2000; 
Lyutikov & Uzdensky 2002). A higher value is justified, however, if the effective internal 
diffusion scale in the current sheet is given by a turbulent mixing length, rather than the 
characteristic Larmor radius of the particles. A turbulent sheet structure arises if the sheet 
filaments into a large number of reconnection centers when field dissipation begins. In the 
model presented here, it is assumed that the global sheet evolution can nevertheless be 
described within an MHD framework, using an effective average value for the dissipation 
rate. 

As the wind approaches the point where dissipation is complete, the plasma temperature 
drops, as can be seen from Figs. 2-5, and the approximation of a relativistic gas fails. One 
effect of this is to invalidate our choice of a fixed ratio of specific heats 7 = 4/3 for the gas 
dynamics. However, more importantly, the dissipation prescriptions we have used are also 
affected. The Harris equilibrium gives a general expression for the drift speed of particles 
in the sheet, valid in both the relativistic and nonrelativistic cases, so that the prescription 
for slow dissipation remains a lower limit. The upper limit given by the limiting the speed 
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of expansion of the sheet to that of sound in an unmagnetized gas can also be implemented 
in the nonrelativistic discussed in Sect. 4.2 and shown in Fig. 5. However, in the 

case of dissipation limited by the growth rate of the tearing mode, the physical basis of our 
estimates is no longer valid once the plasma temperature drops below the particle rest mass. 
An improved treatment, using an expression for the growth rate valid in the transrelativistic 
and nonrelativistic regimes is possible, but our estimates are in any case very rough, and the 
added complications would not appear to be justified. 

However, given that sufficient particles arc present in the wind to satisfy condition (50), 
and neglecting the modification of the tearing mode growth rate for nonrelativistic plasma, 
the prescription for tearing mode limited dissipation automatically lies between the limiting 
dissipation rates and can therefore be considered as a physically realistic possibility. For 
Crab pulsar parameters, L — 1.5 x 10^^ and condition (50) gives fj, < 10^, for Pc — l/'s/S. 
Estimates of the /i parameter for the Crab, based on modelling the X-ray and gamma-ray 
emission (Hibschman & Arons 2001a), exceed this value. In this case the solution obtained 
by LK for slow dissipation loses its validity before dissipation is complete. Nevertheless, 
the range of validity extends from close to the pulsar to well beyond the radius at which a 
termination shock is indicated by observations. X-ray observations (Weisskopf et al 2000) 
place this at a radius of 0.14 pc or r = 2.7 x lO^rL in the equatorial plane. Slow dissipation 
implies steady acceleration of the wind to beyond this radius, without significant dissipation 
of the Poynting flux. 

Recently, however, two different lines of investigation have suggested that the conven- 
tional approach seriously underestimates the rate at which pairs are injected by the Crab 
pulsar into its surrounding Nebula. Gallant et al (2002) have proposed a model that accounts 
for the synchrotron spectrum of the Nebula from the radio to hard gamma-ray bands. In 
particular, the large number of radio-emitting electrons, which could not be explained in the 
model of Kennel & Coroniti (1984), arise naturally from a relatively low speed, high density 
pulsar wind. The estimated injection rate is N± pa 3 x 10^'^ s~^, implying ~ 2 x 10^. 
Lyubarsky & Eichler (2001), who analysed the dynamics of the X-ray jet, also conclude 
that a relatively low value of a at the base of the flow is required. Because of the difficulty 
of coUimating a very low density wind, they suggest that on the equator a < 10^. For a 
transonic wind, this implies /x = 10^. 

Such low values of ^ and the corresponding high pair densities, have important impli- 
cations for dissipation in the pulsar wind. Firstly, independent of the actual mechanism, 
the asymptotic Lorentz factor of the wind after dissipation equals /i. This fits in well with 
the spectral modeUing of Gallant et al (2002), who correlate this value with the break at 
10^^ Hz in the UV spectrum of the Crab Nebula. Secondly, although the radius at which 
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dissipation is complete is independent of ji in the 'slow' dissipation prescription of Coroniti 
(1990), Michel (1994) and LK, this is not true of the faster tearing mode limited case. From 
Table 1 it is evident that for /x = 2 x 10^, dissipation by the tearing mode limited mechanism 
is still not complete within the termination shock. However, a faster mechanism could pro- 
vide the required dissipation in this case. The lower limit on the dissipation radius implied 
by the estimates in Table 1 is 5 x lO^rL and Fig. 5 shows that when this limit is modified 
to take account of the nonrelativistic sound speed, the solution dissipates the magnetic field 
completely at r = 2 x lO^rL, i.e., within the termination shock. Adopting the value ji — 10^, 
or lower, reduces the dissipation radius still further, giving, for the tearing mode limited 
prescription, r = lO^rL. The minimum distance required to dissipate the magnetic field into 
relativistic plasma according to Tabic 1 is r > 1.5 x lO^rL, and a numerical integration using 
the more general limit implied by Eq. (34) gives r > 5 x lO^rL. 

It is interesting to note that the appearance of a dissipative wind may differ between 
these two cases. Kirk et al (2002) derived the criterion 7^ > r/rL for the radiation emitted 
at a given radius in a striped wind to appear pulsed to the observer. Accordingly, for fast 
dissipation, radiation emitted by the hot current sheet at all positions within the flow should 
appear pulsed. On the other hand, the situation for dissipation limited by the tearing mode 
growth rate is more complicated: radiation from close to the maximum radius should not 
appear pulsed. However, nearer to the star the Lorentz factor of the flow is sufficient to 
ensure a pulsed appearance. The magnetic field in the fiow decreases with increasing radius, 
so that the efficiency of conversion of heat into synchrotron radiation should be greater at 
smaller r. Consequently, pulsed radiation can be expected also in the tearing mode limited 
case. 



8. Conclusions 

We model the dissipation of an oscillating, toroidal magnetic field component in a rela- 
tivistic MHD outfiow by considering the radial evolution of the entropy wave associated with 
this pattern. Because of magnetic reconnection, or, more precisely, annihilation, plasma in 
the current sheets separating regions of oppositely directed field is heated. We assign an 
energy density and pressure to the sheets and relate them by the equation of state of a 
relativistic gas. Upon expanding, the sheet performs work, which causes the outfiow to 
accelerate. 

Despite recent progress (Zenitani & Hoshino 2001; Lyutikov & Uzdensky 2002) the 
dissipation rate in a relativistically hot pair plasma embedded in stripes of cold, magnetized 
plasma cannot be written down a priori. Thus, we consider three different prescriptions. 
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each with a specific physical motivation. A lower limit on the dissipation rate is given by 
the prescription used by Coroniti (1990), Michel (1994) and LK, in which the characteristic 
drift speed in the current sheet is equal to the speed of light (to within a factor of order 
unity). We argue that an upper limit on the dissipation rate is set by allowing the sheet 
to grow at its internal sound speed (in the comoving frame). When consistent with these 
limits, a physically plausible estimate of the actual rate is given by assuming that expansion 
of the current sheet is limited by the growth rate of the relativistic coUisionless tearing 
mode. Comparing our prescriptions to a coUisional MHD model, we find that the required 
dissipation is slower than the Bohm value when dissipation is limited by the tearing mode, 
but faster when the sheet grows at its internal sound speed. However, the complex flow 
pattern which will presumably be built up in the dissipation region make it difficult to draw 
conclusions from this analogy. 

Analytic, asymptotic solutions, valid for relativistic plasma temperatures, are found for 
each prescription. These agree with and extend those presented by LK for slow dissipation, 
with 7 oc r^/^, and by Drenkhahn (2002) for fast dissipation, with 7 oc r^/^. For tearing- 
limited reconnection we find 7 oc r^/^^. Numerical integration of the full equation set confirms 
these solutions and demonstrates that they are valid until dissipation is nearly complete, 
although the dissipation rate in the fast prescription slows down once the plasma temperature 
becomes nonrelativistic. In general, transients associated with the initial conditions vanish 
rapidly. However, these transients reveal that for some initial conditions the slow dissipation 
prescription is unphysical in that it fails to create entropy. 

Using the asymptotic solutions, the radius at which the oscillating part of the magnetic 
field is completely dissipated can be estimated in terms of the magnetization parameter /i 
introduced by Michel (1969) and the dimensionless power of the flow L. In the case of the 
Crab, it was shown by LK that slow reconnection cannot dissipate the magnetic energy 
inside the radius at which observations imply that the wind terminates, independent of 
the value of fi. Faster dissipation prescriptions do not change this conclusion, provided 
standard estimates of the rate at which the pulsar injects pairs into the Nebula are adopted 
(Hibschman & Arons 2001a). 

However, our calculations show that a higher injection rate permits faster dissipation. 
Thus, the cr-problem is resolved if dissipation is fast and /i ^ 2 x 10^, corresponding to a pair 
flux of = 3 X 10^'^s^^. This value is consistent with the model of the synchrotron spectrum 
of the Crab nebula presented by Gallant et al (2002). Lyubarsky & Eichler (2001) quote 
an even lower value of ~ 10^ from their considerations of the jet emerging from the Crab 
pulsar. In this case, the cr-problem is also resolved by dissipation proceeding at the tearing- 
mode limited rate. This suggests that most of the energy emitted by the Crab pulsar as 
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Poynting flux is converted into particle-born energy during the supersonic expansion phase 
of the wind. If a fraction of this power is converted into radiation, the dynamics of the flow 
imply that it will appear pulsed to an observer, which lends support to a recent suggestion 
concerning the origin of the optical to gamma-ray pulses from this object (Kirk et al 2002). 

We thank Yves Gallant and Yuri Lyubarsky for stimulating discussions and helpful 
comments on this work. 



A. The relativistic Harris equilibrium 



The relativistic Vlasov equation: 



df df df 



where v is the particle 3- velocity, p = 'ymv the momentum (with 7 = l/y'l — v'^/c^), f the 
Lorentz invariant phase-space density, 

^ = q(^E + ^vAB^ (A2) 

— * — * 

the Lorentz force, q the particle charge and E and B the electric and magnetic flelds, is 
solved by an arbitrary function of the constants of motion. 

In the case of a current sheet in the y-z plane, where quantities in (Al) depend on 
space only via x, the constants of motion are the energy and the y and z components of the 
canonical momentum: 

E = (mV + cV)'^', 

Pz = Pz + -^A,{x), (A3) 

— * — * 

where A is the vector potential. We seek a solution in which the z component of B reverses 
at X = 0. Choosing Az — 0, this means Bz — dAy/dx = at x = 0. To this we can 
add a constant component along y by specifying A^ — zBy, without affecting the following 
analysis. 

Harris (1962) found a solution to this problem with a phase-space density which, at 
X — 0, consists of two Maxwellians, one for each species, drifting in opposite directions along 
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the y axis with zero charge density. The relativistic generahsation found by Hoh (1968) has 
two drifting Jiittner/Synge distributions. In the case with just positrons q = e and electrons 
q = — e, these distributions drift with equal and opposite speeds c^q/e and have equal and 
opposite charge densities. In terms of quantities measured in the rest frame of a species 
(denoted with a bar), the phase-space density of that component at a; = is 



fi 



o± 



/( 



N 



o± 



-exp i-E/e) , 



(A4) 



where E is the particle energy in the rest frame, the temperature and K2 a modified 
Bessel function. The distribution is normalized such that the particle number density is N: 



Employing the Lorentz transformation 

E = T{E^c(3py), 



(A5) 



where the upper sign refers to the positrons, we can write the phase-space density at x — 
in terms of lab. frame quantities: 



/( 



N 



o± 



-exp 



e 



A7rm^cQK2 (mc^/Q) 

Choosing Ay{0) — 0, the solution of the Vlasov equation which matches (A6) at x = is 

N 



(A6) 



/: 



± 



X 



T{ETc(3py^qPAy{x)) 



exp 



/o±exp 



e 



'e(3VAy{x) 

e 



The current density carried by the electrons and positrons is: 

el3TAy{x) 



Jy = exp 
Using the invariance of d^p/E: 



e 



ec' I dV^ (/o+ - /o-) . 



(A7) 



(AS) 



/ 



d'p^fo± 



dy~fo±, 

^2 



= J d'p^r{py±(3E/c)fo±, 
= ± / d^pVPcfoi = ±rpcN, 



(A9) 
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so that 



jy — 2erpcNexp 



e(3rAy{x) 

e 



(AlO) 



Using the Maxwell equation iirjy/c = —dBz/dx we find: 



-87r/3eN^exp 



e 



(All) 



where A^^ = is the number density of each species in the lab. frame. 

With the boundary conditions dAy/dx = ^4^^ = at x = the solution of (All) is 

e 



e(5V 



In [sech^ ■, 



with 



Ad//?, 



Ad 



e 



(A12) 



(A13) 



This is precisely the solution found by Harris (1962), with the substitution T — O/F. The 
distribution function is 



/± = /o±sech^ {x/a) 
and the density and magnetic field are 

= A^^sech^ (x/a) , 

B' = y ^^^tanh(a;/a) 

The number of particles per unit area of the sheet is 



(A14) 



/oo 
dxn'^ix) = 4a at;. 
■oo 



(A15) 



(A16) 



To calculate the elements of the stress-energy tensor 



(A17) 
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where = {E, cp) is the momentum 4- vector, it is convenient to transform to the rest frame 
of each species in turn and use the integrals 



3© + mc 



K2{mcye) 



NR{mc^/e), 
3NQ, 



to find: 



2^1,1 



2^3,3 



r{R + e) 
> © 



2n 



T 



-2,2 



2n' 



± 



r 

r (i? + 9) 



© 
r 



R 

r 



(A18) 
(A19) 

(A20) 
(A21) 
(A22) 



and all remaining elements zero. In the pulsar wind, purely radial motion is considered, with 
the sheet normal along the radius vector. In this case, only the T°'° and T^'^ components are 

relevant and we can identify the effective enthalpy density and prcssiu^c: w = 2n'^r (7? + 0), 
p = 2n^0/r. This leads to an effective ratio of specific heats [see Lyubarsky & Kirk (2001) 
Eq. (A9)]: 



7 



1 + 



© 



r (ri? - mc2) + r2/32e' 



(A23) 



B. The entropy condition Eq. (27) 



Details of the derivation of eqs. (1), (2) and (3) are given in LK. Using the same notation, 
we derive here the condition that the specific entropy of a fluid element always increase, 
Eq. (27). 

The fast phase variable and slow radius variable R, Eq. (A15) of LK are defined by 

dr' 



= Q 



t 



v{r') 



R = er/ri, , 



(Bl) 



where ro is an arbitrary point in the flow. The starting point of the derivation is the entropy 
equation written in terms of the energy density e' and entropy density w' in the co-moving 
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frame, this reads: 



. V \ de' , 



p B?v^ d(j) 



Expanding the dependent variables as in Eq. (A18) of LK and collecting terms of zeroth 
order one finds, using conditions (A19-A22) of LK, that this equation is trivially satisfied, 
because the ohmic dissipative term j' .E' vanishes to zeroth order. To first order, we find an 
expression equivalent to (A27) of LK: 

7oWi de'^ 7o< dvi de'w'^d, 2\ 

- ^|?-I^4(^-oi^o) , (B3) 
vq o(f) A-k^qR oR 

in which the right hand-side is proportional to the ohmic dissipation. To derive the entropy 
equation Eq. (3), this equation is integrated over and exact periodicity is imposed upon 
the expanded dependent variables. However, in order to obtain Eq. (27), we instead impose 
the restriction that the ohmic dissipation be positive: 

7oVi de'o 7oW^ dvi 



Vq dcf) Vq 90 

In order to integrate this condition, we introduce the ratio of specific heats 7 as in Eq. (12) 
and assume this to be constant. Multiplying by the (positive) factor p'~^'^~^^/^^ integrating 
and imposing exact periodicity leads to: 



Noting that for the Harris current sheet 

-1 - . 

Stt / 2 



p'Vn = aH I da:[sech(a;)]'/n , (B6) 



we arrive, after introducing the dimensionless variables used in Sect. 3, at condition (27). 
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C. Expansion speed of the current sheet 

Assuming the sheet is bounded in the comoving frame by a forward and a backward 
moving front at r = r+(r, t) and r_(r, t), respectively, we can write their speeds in the lab. 
frame as 

dr± ^ v±c(3c 

dt l±vpjc ' ^ ' 

where ±/9c arc the speeds of the fronts measured in an inertial frame moving with the speed 
of the wind at (r, t). In terms of the fast phase and slow radius variables, Eq. (Bl), we find 

d0± - (cl/±/(lO iOt/OR) + (Or /OR) ^^2) 



dR -{dr±/dt){dt/d(l)) + {dr/d(f)) ' 

The fraction of one wavelength occupied by the two sheets thus evolves according to 

dR = ndR^"^^-"^-^^ 

'^^^ (C3) 



This is equivalent to Eq. (31), provided r is identified as a 'slow' variable. 
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Table 1. Asymptotic flow properties for different dissipation prescriptions. Given are the 
index q in Eqs. (40-41), the approximate radius rmax at which dissipation is complete, the 
approximate radius rmin at which the solution crosses the fast magnetosonic point and the 
exact (asymptotic) expression for A. The light-cylinder radius is denoted by Tl, A is the 
thickness in units of the pattern half-wavelength, is the total energy carried by the wind 
per unit rest mass (xc^), L is the dimensionless flow luminosity defined in Eq. (29), and (3c 
is the sheet expansion speed in the comoving frame, in units of c. 
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Slow dissipation 
(Eq. 30) 
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Tearing mode limited 
dissipation (Eq. 37) 
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Fast dissipation 
(Eqs. 31 & 33) 
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